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Plastic events in amorphous solids can be much more than just “shear transformation zones” 
when the positional degrees of freedom are coupled non-trivially to other degrees of freedom. Here 
we consider magnetic amorphous solids where mechanical and magnetic degrees of freedom interact, 
leading to rather complex plastic events whose nature must be disentangled. In this paper we 
uncover the anatomy of the various contributions to some typical plastic events. These plastic 
events are seen as Barkhausen Noise or other “serrated noises”. Using theoretical considerations 
we explain the observed statistics of the various contributions to the considered plastic events. The 
richness of contributions and their different characteristics imply that in general the statistics of 
these “serrated noises” cannot be universal, but rather highly dependent on the state of the system 
and on its microscopic interactions. 


I. INTRODUCTION 

Modeling the mechanical properties of amorphous 
solids is an active subject of current research, requiring 
detailed understanding of the many-body processes that 
occur in such system when subjected to external strains. 
External strains can be mechanical, magnetic or electric, 
depending on the properties of the amorphous solid in 
question. The responses of amorphous solids to such ex¬ 
ternal strains is usually not smooth, giving rise to ’’ser¬ 
rated” plots of stress vs. strain, energy vs. strain, mag¬ 
netization vs. external magnetic held etc. A lot of effort 
was spent on characterizing the probability distribution 
functions of such serrated responses. In the context of 
magnetic jumps this is referred to as Barkhausen Noise 
[iHg but other serrated noises were studied as well Q- 
In a number of cases strong claims of universality were 
made. 

Recently we have analyzed in some detail model amor¬ 
phous solids in which there is a signihcant coupling be¬ 
tween mechanical and magnetic properties |8l-lll| . Doing 
so we realized that the characterization of the physics 
of plastic events can be quite demanding; there is more 
in these events than what meets the eye at hrst impres¬ 
sion. The aim of this paper is to highlight the some¬ 
what complex anatomy of plastic events in such systems, 
with a word of caution to researchers in the field that 
similar complexity may arise in other systems as well, 
and reasonable modeling should take this into account. 
In particular we will conclude below that in general one 
should not expect universal probability distribution func¬ 
tions since the statistics of the serrated responses depend 
on many details of the microscopic interactions and on 
the state of the system. 

Deferring all details to the next section, we motivate 
the present paper by showing in Fig. [T]a scatter plot of 
the values of energy drops AU during plastic events when 
the system is strained by an external magnetic field. The 
scatter plot is shown as a function of the magnetization 
jump Am that occurs simultaneously with the energy 
drop. First, one sees that for a given Am one has a wide 
distribution of AU values. Second, these values of the 
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FIG. 1. Scatter plot of value of the energy drops AU as a 
function of the simultaneous changes Am in the magnetiza¬ 
tion. These occur when the magnetic field is ramped up and 
down to form the hysteresis loop of the Barkhausen Noise. 


energy drops fall in different groups, with a strange in¬ 
tense line and two triangular groups that are only partly 
overlapping. Understanding such scatter plots and their 
implications on the physics of the solid is what we mean 
by “the anatomy of plastic events”. A theory of plasticity 
in such system should include also the understanding of 
the statistical distribution of such events. The density of 
points in every little box of size dAUdAm inFig. [T]is pro¬ 
portional to the joint probability P{AU, Am)dAUdAm. 
Understanding how this probability distribution function 
is determined by the different physical process requires 
a theory of the anatomy of plastic events. The aim of 
this paper is to provide such a theory for one particular 
model of amorphous magnetic glass. Other models will 
require a similar approach. 

In Sect. m we present the model that was introduced 
recently and analyzed for some of its aspects in Refs. [1- 
[TTI| . The next section Hill deals with the notion of plastic 
events, both under external mechanical strain and un¬ 
der external magnetic field. The anatomy of the plastic 
events that occur under magnetic straining is studied in 
section IlYl in which we provide the numerical analysis, 
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a theory, and a comparison between the two. In Sect, 
m we discuss briefly the anatomy of plastic events under 
mechanical straining, including the implications to mag¬ 
netostriction. The last section offers a summary and a 
discussion. 

II. A MODEL OF A MAGNETIC AMORPHOUS 

SOLID WITH STRONG LOCAL ANISOTROPY 

The model Hamiltonian was introduced in Q and an¬ 
alyzed further in [Ml- The magnetic part of the model 
Is in the spirit of the Harris, Plischke and Zuckerman 
(HPZ) Hamiltonian [l^ but with a number of important 
modifications. These modification were made to bring 
the model closer to the physics of amorphous magnetic 
solids 0. The first major difference is that the particles 
In the present case are not pinned to a lattice. We write 
the Hamiltonian as 

C/({n}, {Sj}) = C/i„ech({n}) + tfmag({u}, {‘S'J) , (1) 

where {ri}fLi are the 2-dimensional positions of N par¬ 
ticles in an area and Si are spin variables. The 
mechanical part Umech is chosen to represent a glassy 
material with a binary mixture of 65% particles A and 
35% particles B, with Lennard-Jones potentials having 
a minimum at positions cTyi^ = 1.17557, cfab = 1-0 and 
<^BB = 0.618034 for the corresponding interacting par¬ 
ticles M- These values are chosen to guarantee good 
glass formation and avoidance of crystallization. The en¬ 
ergy parameters chosen are caa = ess = 0.5 eAB = 1-0, 
In units for which the Boltzmann constant equals unity. 
All the potentials are truncated at distance 2.5cr with 
two continuous derivatives. JVa particles A carry spins 
Si; the Nb B particles are not magnetic. Of course 
Na + -IVb = N. We choose the spins Si to be classi¬ 
cal xy spins; the orientation of each spin is then given by 
an angle (pi with respect to the direction of the external 
magnetic field which is along the x axis. 

The magnetic part of the potential energy takes the 
form Q: 

<i]> 

- ^ Ki cos^ {pi - 6»i({ri})) - y,AB ^ cos {pi) . (2) 

i i 

Here = \ri — rj\ and the sums are only over the A 
particles that carry spins. Notice that in the present 
model the exchange parameter J{rij) Is a function of a 
changing inter-particle position (either due to affine mo¬ 
tions induced by an external strain or an external mag¬ 
netic field or due to non-afRne particle displacements, 
and see below). Thus randomness in the exchange in¬ 
teraction is coming from the random positions {r^}, 
whereas the function J(Uj) is not random. We choose 
the monotonically decreasing form J{x) = Jof{x) where 
f{x) = exp(—x^/0.28) + Hq + H 2 x‘^ + H^x'^ with Hq = 


-5.51 X 10"® ,772 = 1-68 x 10“® ,774 = -1-29 x 10"®. 
This choice cuts off J{x) at x = 2.5 with two smooth 
derivatives. Note that we need to have at least two 
smooth derivatives in order to compute the Hessian ma¬ 
trix below. Finally, in our case Jq = 3. 

Another major difference with the HPZ model Is that 
In the present case the local axis of anisotropy 9i Is not se¬ 
lected randomly, but Is determined by the local structure. 
Recall that In a crystalline solid the easy axis is deter¬ 
mined by the symmetries of the lattice. In an amorphous 
solid the arrangement of particles changes from one posi¬ 
tion to the other, and we need to find the local easy axis 
by taking this local structure Into account. Define the 
matrix 

Tf ^ Y. k4 / E ) • (3) 

3 3 

Note that we sum over all the particles that are within the 
range of J {rij ); this Is sufficient to take Into account the 
arrangement of the local neighborhood of the ith particle. 
The matrix Ti has two eigenvalues in 2-dimensions that 
we denote as and Kip, Kip > Kip. The eigenvector 
that belongs to the larger eigenvalue Kip Is denoted by n. 
The easy axis of anisotropy Is given by 9i = sin“^(|ny|). 
Finally the coefficient Ki which now changes from parti¬ 
cle to particle Is defined as 

77i = J(rij)]^(Kip-Kip)^ , C = Ko/Jo(Jab- (4) 

3 

The parameter Kq determines the strength of this ran¬ 
dom local anisotropy term compared to other terms in 
the Hamiltonian. For most of the data shown below we 
chose Kq = 5.0. The form given by Eq. ensures that 
for an Isotropic distribution of particles Ki = 0. Due to 
the glassy random nature of our material the direction 9i 
is random. In fact we will assume below (as can be eas¬ 
ily tested in the numerical simulations) that the angles 
9i are distributed randomly In the interval [—7r,7r]. It Is 
important to note that ramping the magnetic field does 
NOT change this flat distribution and we will assert that 
the probability distribution P{9i) can be simply taken as 

P{9,)d9, = ? . (5) 

ZTT 

We have checked In the numerical simulations that Eq. 
([5]) is valid to a high approximation at all values of B. 
The last term In Eq. ([J) is the interaction with the exter¬ 
nal field B. We have chosen /rpTl in the range [-0.08,0.08]. 
At the two extreme values all the spins are aligned along 
the direction of B. 

In passing we should comment on the chosen param¬ 
eters in the model. Our guiding line was to choose pa¬ 
rameters such that the magnetostriction coefficient Is of 
the order of what Is known in laboratory materials. De¬ 
creasing Kq results In much smaller magnetostriction co¬ 
efficients and vice versa. We did not aim at modeling a 
particular material, and our Interest here, as before. Is 


ill the generic properties of amorphous solids with strong 
local anisotropy. 
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III. WHAT IS PLASTICITY? 

It is well known that plasticity in crystalline solids is 
carried by defects, like dislocations, whose glide under ex¬ 
ternal strains is dissipative, leading to energy loss. What 
are the mechanisms of energy loss in amorphous solids is 
less well known, although research in the last two decades 
has shed considerable light on the fundamental physics 
of plasticity in amorphous solids In the present 

system we can have two distinct external agents that can 
strain the system, i.e. mechanical strain and magnetic 
field. Such external strain can be studied in systems hav¬ 
ing finite or zero temperature. Since we are interested in 
the anatomy of plastic events we opt for the latter, tem¬ 
perature fluctuations tend to mask the clear cut plastic 
events that are recognized at T = 0. To keep the sys¬ 
tem at T = 0 we must also ramp the external strain or 
the magnetic field quasistatically to allow the system to 
remain in mechanical equilibrium at all times, without 
heating effects. In such conditions it is completely clear 
what are the instabilities that are responsible to plastic 
events. 

The response of our system to external strain, be it 
mechanical or magnetic, is reversible and smooth as long 
as the system is mechanically and magnetically stable. 
This is the case as long as the Hessian matrix T-L has 
only positive eigenvalues. In the present case H takes on 
the form Q: 


n = 


d^U 

dvidrj 

d^U 

d<pidri 


d^U \ 

dvidcpj \ 

a^u I 

dcpidcpj / 


( 6 ) 


The system loses stability when at least one of the eigen¬ 
values of "H goes to zero. When this happens, there ap¬ 
pears an instability that results in a discontinues change 
in stress, in energy and in magnetization. In Fig. [2] we 
show a typical blown up section of the energy, magne¬ 
tization and stress curves as a function of B. We see 
that the discontinuities appear simultaneously in all the 
three quantities at the same values of B. These are ir¬ 
reversible plastic events that take the system from one 
minimum in the energy landscape through a saddle-node 
bifurcation to another minimum in the energy landscape 
where again all the eigenvalues of H are positive. In Ref. 
Q we derived an exact equation for the dependence of 
any eigenvalue \k on B for a fixed external strain, which 
reads: 


I _ (b) _ [bkM + 

dB ^ \t 

t 


The precise definition of all the coefficients is given ex¬ 
plicitly in Ref. Q . Generically, when one eigenvalue, say 



FIG. 2. The magnetization, the energy per particle and the 
stress component as a function of external magnetic field 
B. This figure demonstrates that all three quantities have dis¬ 
continuities at the same value of B where the system under¬ 
goes a plas'' ' . ...1 . ■ n., . TT .... 

matrix n ] 
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FIG. 3. The logarithm of the eigenvalue Ap that hits zero at 
Bp as a function of the logarithm of Bp — B. The slope has 
a value of 1/2. 


Ap approaches zero, all the other terms in Eq. o remain 
bounded, leading to the approximate equation 

9Ap, Const. , , 

® 

In such generic situations the eigenvalue is expected to 
vanish following a square-root singularity, Ap ^ {Bp — 
where Bp is the value of the external magnetic 
field where the eigenvalue vanishes. The reader should 
be aware of the fact that at some special values of B it 
may happen that the coefficient Const in Eq. [8] vanishes 
at the instability leading to a an exponent different from 
1/2 i. This non generic feature hardly changes the con¬ 
siderations of the present paper. In Fig. [3] we show a 
typical dependence of the eigenvalue Ap on B, where the 
square-root singularity is apparent. It is also interest¬ 
ing to examine what happens to the eigenfunctions 
which are associated with the eigenvalues Xk as the in- 
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FIG. 4. The projection of the eigenfunction 4'^ associated 
with the eigenvalue Xp shown in Fig. [3] projected on the par¬ 
ticles positions and on the spins in the upper and lower pan¬ 
els respectively. The upper panel shows a typical non-afhne 
displacement field associated with a plastic event, having the 
quadrupolar structure of an Eshelby solution. The lower panel 
shows that the same event is associated with a co-local flip of 
spins, leading to the change Sm of the Barkhausen Noise. 


stability is approached. The answer is that all the eigen¬ 
functions of H are delocalized far from the instability, 
but the one eigenfunction associated with Ap —>■ 0 
gets localized on n N particles. A typical projection 
of close to the instability on the particles positions 
and on the spins is shown in the two panels of Fig. 0] We 
see that the non-afhne movement of the particles is very 
similar to the standard “Eshelby like” quadrupolar event 
that is so typical to amorphous solids. The projection on 
the spin degrees of freedom shows that a patch of spins 
had changed its orientation (magnetic hip of a domain). 
Note that the patch is compact, without any fractal or 
other esoteric characteristics that were associated with 
Barkhausen Noise in the past. This is the nature of the 
event that is associated with the Barkhausen Noise in 
our case. The reader should be aware however of the fact 
that the addition of long range dipole-dipole interactions 
can change this qualitatively, leading to elongated mag¬ 


FIG. 5. Scatter plot of value of the stress drops Aa as a func¬ 
tion of the simultaneous changes Am in the magnetization. 
These occur when the magnetic field is ramped up and down 
to form the hysteresis loop of the Barkhausen Noise. 

netic domains and a different mechanism of Barkhausen 
Noise due to the movement of domain boundaries 


IV. THE ANATOMY OF PLASTICITY 

The interesting physics of plasticity in this model stems 
from the fact that the Hessian matrix (jS]) couples the po¬ 
sitional to the magnetic degrees of freedom. Thus a plas¬ 
tic drop in stress and energy will be usually coupled also 
to a change in the magnetization. Whether one strains 
the system with a mechanical strain or a magnetic field, 
the plastic drops will be composite processes in which all 
the degrees of freedom contribute to the non-afhne re¬ 
sponse. In this section we focus on these events that are 
triggered by the magnetic held as the straining agent. To 
expose the anatomy of the plastic events we re-plot the 
data in Fig. [2] in a different way, i.e. as a scatter plot 
of the drops in energy or in stress as a function of the 
magnetization change Am. The hrst was shown in Fig. [1] 
and the second is shown here as Fig. [3 As before, we see 
that also the values of the stress drops organize into two 
distinct groups that are however not non-overlapping. 


A. Detailed analysis of the energy drops 

Our hrst task is to rationalize the distributions that 
appears in hgures like Figs. [Hand[H Focussing as an ex¬ 
ample on the energy drops, we return to the Hamiltonian 
and hnd out which of the terms is responsible to which 
group of energy drop values in these hgures. This sepa¬ 
ration is demonstrated in Fig. [51 where the energy drop 
is assigned to four different contributions to the Hamil¬ 
tonian, i.e. the Lennard-Jones positional degrees of free¬ 
dom, the exchange interaction, the anisotropy energy and 
hnally the interaction with the magnetic held. Obviously 
the combination of the scatter plots in Fig. [5] will lead to 
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FIG. 6. The four distinct contribution to energy drops in plas¬ 
tic events. The combination of all these scatter plots should 
yield the data in Fig. [T] 


what was shown as Fig. [TJ The same decomposition can 
be done for the stress drops but for the sake of brevity 
we focus here on understanding the results shown in Fig. 
[51 a similar analysis for the stress drops is implied. 

To understand what we see we will invoke the re¬ 
sult of the previous work [TTI | in which the distribu¬ 
tion P{Am) was measured over three orders of mag¬ 
nitude 10“^'® < Am < 10“ ®. These magnetization 
changes involved flips of magnetic domains of between 
10° < An < 10° particles. This probability distribution 
function (PDF) was found to be well fitted in this regime 
by a form 


we can assume a Gaussian form for the conditional prob¬ 
ability 


P2(AC7|Am) = exp-(At/-(At/))V(2g^). 

v27r cr° 


(13) 

where both the average energy drop (AC/) = (AC/) (Am) 
and the variance of the energy drops cr^ = tT(Am)° are 
functions of Am. 

We can now introduce two exponents and C 2 by 


(AC/) (Am) = KiAm^^ 

cr^(Am) = K 2 Am‘’^. (14) 


One of the purposes of these notes is to estimate these 
exponents Ci and C 2 . 

Now the total energy drop can be separated into its 
individual mechanical and magnetic contributions 

AU = AC/’"“'" -b AC/™®® 

= AC/’"“'^-b AC/“-b AC/““'’-b AC/^ (15) 

where AC/™“'" is the mechanical contribution to the 
energy drop, while the spin contribution can be sepa¬ 
rated into exchange AC/®“, anisotropic AC/“"*® and mag¬ 
netic field AU^ contributions, and depending on the spin 
Hamiltonian for the metallic glass (which in our case is 
Eq. ([2])) may lead to different exponents. 

Let us therefore analyse the consequences of Eq. (ED 
carefully. First we note that 


(AC/)(Am) = ^(u,), (16) 

2=1 


P,A„) = , (9) 

where the exponential decay rate A is analytically com¬ 
puted and the function /(Am) is evaluated explicitly. 

To make the connection to our present data we need to 
discuss the conditional distribution P 2 {AU\Am) in terms 
of which the joint distribution 

P(AH, Am) = P(AC/| Am)P(Am) (10) 


and thus if a non zero {ui) = (u) exists we would find that 
(AC/) (Am) ^ Am or /i = 1. On the other hand if there 
exists a contribution consisting of n random variables 
with zero mean, then we might expect that contribution 
to scale like ~ Am°/^ or Ci = 1/2. The important point 
to note is that by measuring the physics underlying 
the spin flips can be found. 

To estimate the variance of the fluctuations cr^ = 
(AC/^) — (AC/)° we write using Eq. (fTT|) 


can be written. 

We shall start with the general form for the energy drop 

n 

AU = Y,u,, ( 11 ) 

i=l 

where Ui is the energy change (both mechanical and mag¬ 
netic) associated with the slip of the fth spin in the flip¬ 
ping domain and n is the number of flipping spins. We 
assume that each spin flip contributes to Am and thus 

n ^ CAm , (12) 

with some unknown constant C. Thus AC/ is a sum of 
n random variables, and using the central limit theorem 


cr(Am)2 = - {ui){uj)]. (17) 

i=l j=l 

Thus while (AC/) only depends on the additive contribu¬ 
tion of n random variables, the variance also depends on 
the correlation of different spins within the flipped do¬ 
main. To see the consequence of these correlations, let 
us consider first the case where the spins i and j are un¬ 
correlated. In that case a{Am)‘^ = n[{u^) — (u)^] ~ Am 
or ((2 = 1. On the other hand, in the limit of strong cor¬ 
relations between spins in the flipped domain [{uiUj) — 
{ui){uj)] yC 0 and as a consequence a{Am)‘^ ^ Am? or 
C 2 = 2. It is also possible that the clean scaling described 
above may not exist but rather several mechanisms are 
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0.5 

1.0 

0.5 

1.0 

C 2 

0.75 

1.0 

1.0 

1.0 


TABLE I. The best fit values of and ("2 for different energy 
terms. Note that the exponent of 0.75 is not explained theo¬ 
retically but the data for this particular contribution is small, 
close to the noise level. 


mixed. Only simulations and data analysis can answer 
these questions. 

At this point we want to use the results of Ref. [ll| 
in which an analytic form for the probability to see a 
magnetic jump of Am was proposed. To this aim we 
will attempt to find analytic approximants to the pdf 
P{AU\Am) in the form of four distinct contributions in 
agreement with Fig. [51 

P(AC7|Am) = PLj{AU\Am) + PExiAU\Am) 
+PAniAU\Am) + PB.siAU\Am) . (18) 



Analyzing the data shown in Fig. |6]we could approximate 
each of the four conditional pdfs in the form 

Pi(A[/|Am) Ri e kTIaToTT ^ 

The best fits for the exponents Ci and (2 and for the 
constant <t are provided in table HI 


The quality of the fits are demonstrated in Figs. [7] and 
IHl The exponents agree with our expectations of being 
1/2 or 1, except for the variance of the Lennard-Jones 
contribution in Fig. jS] We note however that the ampli¬ 
tude of the Lennard-Jones AU is two orders of magnitude 
less than the other contributions and therefore we we are 
close to the noise level and cannot trust this particular 
measurement. The smallness of this contribution arises 
from the fact that the straining here is done magneti¬ 
cally and the position of the particles do not change that 
much. 

Finally we write P{AU) in the form 


P{AU) 


= / dAm 


PLj{AU\Am) + PEx{AU\Am) 



FIG. 8. The variance in the energy changes as a function of 
Am for the four different contributions that play a role in our 
system. 


+PAn{AU\Am) + PB.s{^U\Am) P{Am) . (20) 


V. STRAINING MECHANICALLY 


Now we use the analytic form of P(Am) from Ref. [TT| 
and compute the integral (I20[) numerically. The compar¬ 
ison of this reconstruction of the pdf of AC/ to its direct 
numerical calculation is shown in Fig. |9l 

The conclusion of this exercise is that providing the 
anatomical details of the plastic events can help in un¬ 
derstanding the statistics of energy or stress drops. We 
do not repeat in this paper the exercise for the stress 
drops since it follows verbatim the same steps. 


A similar richness in the anatomy of plastic events 
is found when the system is strained mechanically M- 
Even though we strain mechanically the coupling be¬ 
tween positional and spin degrees of freedom results again 
in having a change in magnetization together with drops 
in energy and in stress. In Figs. llOl and llll we show a typ¬ 
ical plot of energy, stress and magnetization vs external 
strain, for zero magnetic field and for a finite magnetic 




















FIG. 9. Comparison of the direct calculation of P{AU) in 
plastic events due to straining with the magnetic field to the 
reconstruction of the same quantity from the anatomical dis- 
S6Ction tViiQ miflntitAr anrl tViP inHpripnHpnf Vnmsflprl d’p nf 

P{Am) Q 
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FIG. 10. Typical dependence of the energy, stress, and mag¬ 
netization for zero extenral magnetic field. Note that all the 
plastic events occur simultaneously for all quantities. In the 
present case the magnetization is fluctuating up and down 
around a zero mean value. 


field. Note that in the first case the total magnetization 
remains zero on the average, with the flips in magneti¬ 
zation Am being negative or positive with equal proba¬ 
bility. For finite magnetic field magnetization is accumu¬ 
lated in the direction of the magnetic field at each plastic 
event. As before with ramping the magnetic field we see 
that also with mechanical strain the plastic events cou¬ 
ple mechanical and magnetic degrees of freedom. The 
drops in energy are occurring as a result of plastic in¬ 
stabilities at the same values of 7 as the drops in stress 
and magnetization. The mechanism is the same, i.e. an 
eigenvalue of the Hessian matrix hits zero punctuating 
the smooth curves of energy, stress or magnetization with 
sharp drops of irreversible events. 



FIG. 11. Typical dependence of the energy, stress, and mag¬ 
netization for external magnetic field B — 0.01. In the present 
case the plastic events cause the magnetization to increase; we 
called this phenomenon “plasticity induced magnetization” 

[ii. 


To understand the serrated response curve of the en¬ 
ergy or the stress one needs again to search for the 
anatomy of the events, displaying carefully the contri¬ 
bution of each physical mechanism for either energy or 
stress drop. Since we did not compute independently 
the ’Barkhausen Noise’ P{Am) in this case we do not 
repeat the exercise for the case of mechanical straining. 
We stress however that any interested researcher must 
pay attention to the rich physics that is underlying the 
serrated ’’noisy” character of the data shown in Figs. fTOl 
andini 


VI. SUMMARY AND CONCLUSIONS 


The main conclusion of this paper is that characteriz¬ 
ing and understanding the statistics of serrated noise is 
not necessarily the same. Even if we can plot the pdf’s of 
energy drops or of magnetic jumps and measure the ex¬ 
ponent that is associated with their log-log plot, it does 
not mean that we uncovered the intricate physics that 
underlies the phenomenon. We have seen here that even 
the simplest coupling between mechanical and magnetic 
degrees of freedom results in a multitude of contributions 
to the energy changes upon plastic events. Each contri¬ 
bution comes with its own statistics, its own exponent 
and its own amplitude. Of course, once we have the full 
information of all the contributions we can reconstruct 
the pdf of any wanted quantity, cf. Fig. O The full infor¬ 
mation is however not always available in experimental 
systems. Thus great care is called for interpreting the 
observed statistics of serrated noises. In particular we 
should stress that changing conditions (like zero or non¬ 
zero magnetic field in Figs. (TU] and ITT]) may change the 
statistics of the serrated noise. The amplitude of the var- 
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